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ABSTR.ACT 

In the present work an attempt has been made to 
investigate the free vibration response of hybrid 
laminated composite plates made of Graphite /Epoxy and 
Kevlar/Epoxy laminae* Higher order ^ear deformation 
theory has been employed to take into account the trans- 
verse shear and rotatory inertia effects on the natural 
frequencies. Rayleigh-Ritz energy technique has been 
employed to solve vibration problem. Convergence studies 
have been carried out by considering number of terras in 
the displacements and rotations. 

The validity of the present formulation has been 
checked by comparing the non-dimensional frequency parameter 
values with those obtained by other investigators for 
isotropic and regular composite laminates. Symmetric as 
well as antisymmetric angleply and crossply plates with 
simply supported edge conditions have been investigated. 

The effect of various parameters such as nuBfcer of layers, 
fibre orientation, aspect ratio, length to thickness ratio 
etc. have been studied and the results are presented for 
regular as well as hybrid composite plates. 





The Eigen-value problem has been solved using a 
standard subroutine available in DEC-1090 system. It 
has been found that the transverse shear effects cannot 
be neglected for plates having length to thickness 
ratio less than or equal to 20. The frequency parameter 
increases with length to thickness ratio thus proving 
that transverse shear effects can be neglected for thin 
plates without any appreciable loss in the accuracy of 
the results . The frequency parameter value increases * 
with aspect ratio and number of layers. This attains an 
asymptotic value as the number of layers i®.liisre<a»©d« 
Frequency parameter has highest value for fibre orientations 
of + 45^ since the stiffness of the plate is maximum for 
these fibre orientations. It has been found tiiat hybrid 
laminates give intermediate properties as compared to 
those of single material laminates. It has been observed 
that hybrid laminates give values almost nearer to those of 
the Grafrfiite /Epoxy laminates for aspect ratios up to 1. 

For aspect ratios upto 3, the frequency parameter reaches 
asymptotic value beyond length to thickness ratio of 
30, while for higher aspect ratios, transverse shear effects 
are found to be predominant for length to thickness ratios 
higher than 30 also. 
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CHAPTER 1 

INTRODUCTION AND LITERATURE REVIEW 
1 *1 Introduction 

Over the last two to three decades the use of 

composite materials consisting of either continuous or 

discontinuous fibres dispersed in a suitable metallic or 
matrix 

nonmetallic Material has received wider attention. Due 
to the inherent tailoring of properties of these mater|.als 
a number of unique design features can be utilised, inclu- 
ding such potential aspects as reduction in weight, 
improved performance, increased service life, reduced 
system maintenance and active controlled configuration of 
structures. In aircrafts controlling the deflected shape 
of the wing during flight is of paramount importance. 

Structural elements such as plates have found wide 

spread application in aeronautical, automobile, civil, 

marine, mechanical and other fields of engineering. 

They are subjected to various loading environment. It 

becomes imperative to have a thorouc^ knowledge of liie 

one 

structural response, so as to enable* to achieve satis- 
factory design of these elements. 

Laminated plates have been extensively used in 
the aerospace engineering industry. In laminated composite 
plate construction, higSi strength, low density fibres are 
dispersed at desired orientations in a matrix of resin* 
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Such a construction results in heterogeneous, anisotropy 
and the behaviour of these plates becomes highly compli- 
cated. To exploit the advantages of anisotropy, one 
should thoroughly investigate, stmctaral response d«iie 
to transverse inplane and dynamic loads. 

Extensive research has been carried out over the 
past two decades and enough information has been documented 
on structural behaviour of regular composite plates. It is 
well known that structural elements made of two or more 
different materials known as h^rids will have improved 
properties over the former. This is the easiest way 
for obtaining the intermediate properties. Damping 
properties of the structure can also be improved by a 
proper choice of constituent materials. 

The objective of a designer is to control iiie 
vibration ^en it is not desirable. If not controlled, 
this may eventually lead to a failure due to fatigue. 

A study of the available open literature reveals 
that II very few investigations have been made on ihe 
dynamic response of hybrid laminated composite plates. 

In most of -die previous investigations either the 
transverse shear effect was neglected, or the first order 
shear deformation theory has been employed. Some 
investiga^MFlJ have used a shear correction factor to take 
into account the transverse shear effects. 
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In the present work an attempt has been made to 
obtain the free vibration characteristics of rectangular 
hybrid laminated composite plates taking into account 
shear deformation and rotatory inertia effects o A 
higher order idieory proposed by Reddy*^ has been employed 
to obtain natural frequencies and mode shapes, using 
Rayleigh-Ritz technique* 

1 *2 Literature Review 

The available literature which are relevant to 
the present work can be classified into the following 
categories: 

1. Classical thin plate theory* 

2* First order shear deformation theory, 

3. Second order shear deformation theory* 

4. Higher order shear deformation theory* 

In what follows, an attempt has been made to explain how 
each of these theories have been used in the past for the 
analysis of isotropic as well as anisotropic plates. 

1*2*1 Classical Thin Plate ITieory 

1 

It was Lekhnitskii , who applied the Kirchoff*s 
hypothesis for the analysis of synnaetric laminates. He 
published a book vrtierein he has en^loyed -the single layer 
anisotropic plate -theory based on classical plate "theory 
for systematically analysing the bending, stability and 



vibration characteristics of specially orthotropic 
laminates. A generally orthotropic plate has been 

2 

analysed using the Kantrovich first iteration method • 

3 

Reissner and Stavsky improved this analysis by 
taking into account bending-extensional coupling in 
the case of unsymmetrically laminated plates. 

4 

V/hitney and Leissa used this theory for ttie dynamic 
analysis of composite plates. They obtained a closed form 
solution for the linearised equations for the cases of 
bending, vibration and buckling of some important classes 
of laminates for which coupling between bending and 
stretching is present. They found that coupling increased 
deflections vrfiile decreased the fundamental frequency 
and buckling loads. They concluded that the amount of 
reduction of the effective stiffness of the plates depends 

(c 

on the degree of anisotropy of the individual levels and 

on the total number of plies in the laminate. 

5 

Durvasula and Srinivaasan using Huff ington and 
Hoppmann’s equations obtained solutions for exact 
frequency equations for orthotropic plates with one pair 
of opposite edges simply supported vhile the other pair 
of edges supported in any combination. They employed 
Rayleigh-Ritz technique for the study of vibration of 
uniform rectangular orthotropic plates with different 
boundary conditions using beam characteristic functions. 
They also found out the buckling parameters for a plate 
under uniform compression using the corresponding 
frequercy parameters. 


Lin studied the free vibration characteristics 
of unsymmetrically laminated plates using classical 
Kirchoff ’s assumptions. For more general boundary 
conditions Bolotin’s assyraptotic method has been 
employed and natural frequencies and natural modes of 
isotropic homogeneous plates and orthotropic rectangular 
plates with simply supported and clamped edges have been 
obtained. 

8 

Bert and Mayberry used the classical thin lamination 
theory for the free vibration analysis of rectangular 
unsymmetrically laminated anisotropic plates with clamped 
edges. They employed the Rayleigh-Ritz technique and 
compared the fundamental frequency obtained using this 
theory with the experimental results. It was found that 
results were comparable with the experimental results 
only for the fundamental mode. 

f 

Patra and Iyengar ' studied the free vibration 
response of laminated rectangular plates with and without 
cut outs after suitably modifying the thin plate bending 
element used for isotropic plates. They observed that 
frequency changes with fibre orientation and that the 
mode shapes continuously change with fibre orientation. 

1 # 11 ) 

IMaretiya and Iyengar * have studied the free 
vibration characteristics of hybrid laminated composite 
plates. They have considered rectangular as well as 
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skew plates vdth one pair of opposite edges simply 
supported and the other pair of edges clamped. They 
employed the Galerkin-Vaslov technique A linear 
displacement field has been considered. They have 
obtained optimal layup details with constraints on 
maximum deflection and fundamental frequency. In 
this investigation it was shown that the h^rid 
laminates give intermediate behaviour between that of 
single material composites \rfiich may be desirable from 
other considerations. 

1 .2 .2 First Order Shear Deformation Theories 

12 

It was Reissner , vtio considered the effect of 

shear deformation on the bending of elastic plates. 

13 

M.ndlin based on similar theory proposed a displacement 
field to take into account the effects of transverse shear 
and rotatory inertia and obtained solutions for the 
flexural vibration of isotropic elastic plates^ ilith 
the following displacement field* 


u 

tl 

X 

X 

y»t) 


= Z «y (X, 

y»t) 

¥ 

= if (x, y. 

t) 


Some investigators have sttidied Hie problem employing 
a shear corroition factor. 

Yang et al proposed another form of shear 
deformation theory. They obtained two dimensional theory 
of motion of heterogeneous plates from three dimensional theory 
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Of elasticity. Transverse shear and rotatory inertia 
was included in the general theory. They proposed the 
following form of displacement field, 

u = UQ(x,y,t) + (x,y,t) 

V = v^(x,y,t) + Z’i' y (x,y,t ) 
w = w(x,y, t) . 

They showed that the frequency curves obtained by using 

classical lamination theory deviated very much from the 

ones obtained using this theory for anisotropic plates. 

This -Uieory is commonly known as YNS theory. Whitney 
IS 

and Pagano > extended this theory for the static and 
dynamic analysis of laminated plates consisting of 
arbitrary number of bonded anisotropic layers. 

17 

Levinson and Muriby - used the equilibrium 

equations of the first order shear deformation -Uieory 

1S IS 

proposed by Whitney and Pagano . Later on Reddy 
showed that these equations were variationally inconsistent 
for the displacement field used, with ibose derived from 
the principle of virtual displacements. 

Bert and Gh®i using the YMS theory obtained 
closed form solution for the vibration of antisymmetric 
angle ply laminated rectangular plates taking into account 
the effect of shear deformation, parametric effects such 
as aspect ratio, length to thickness ratio, nuiBber of layers 
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and fibr® orientatipn •simply supported boundary conditions 
on all edges are considered. It was found that aspect 
ratio effect was more pronounced for thick plates than 
for thin plates and that error in classical plate theory 
as compared to shear deformation theory increases with an 
increase in the number of either longitudinal or transverse 
waves. 

20,21 

Patra and Iyengar used a displacement field very 
much similar to the one used by Mindlin for the free 
vibration study of composite plates by finite element method. 
The displacement field throucfrj out the plate is as belows 

u(x,y,z) = 

v(x,y, 2 ) = -Z (x,y) 

w(x,y, 2 ) = 

where 0„ and 0^, are the rotations of the normals and * 

X y 0 

is the transverse displacement of the raid plane. They 
used a 8-noded, 3 degree of freedom finite elements for 
this, study. 

22 

Reddy studied the free vibration of antisymmetric 

angle ply laminates including transverse shear deformation 

19 

effects using Finite Element Method. Like Bert and C3ien 
he also used ttie YNS theory. Non-dimensional frequency 
parameters for isotropic plates were attained and compared 
with the solutions obtained from 3-D linear elasticity. 
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Results obtained for antisymmetric plies were compared 

1«9 

with the closed form solutions of Bert and Chen . 

It was found that for a given aspect ratio finite 
element results approached closed form results as side 
to thickness ratio decreased. A(2x2)mesh of quadratic 
element was used for quarter plate in all these studies. 

23 

Pryor and Barker ' have used a finite element 
model including transverse shear deformation effects for 
arbitrary laminated plates. They have employed a 28 
degree freedom discrete element* The displacement 
formulation used corresponds to first order shear deformation 
theory. They have used this only for the static analysis 
of laminated plates and have ^own the effect of transverse 
^ear deformation on deflelstion and stresses in comparison 
to the classical lamination theory. 

24 

Joshi and Iyengar employing finite element 
technique have carried out the static and dynamic analysis 
of h'i^rid laminated composite plates. The shear 
flexibility is included in the finite element modelling. 

A nine noded isoparametric element has beei employed 
using Whitney and Pagano’s laminate theory wherein 

transverse shear and rotatory inertia have been taken 

t3 

into account as in Mindlin*s theory for isotropic plates* 



The non-dimensional frequency and deflection parameters 

thick and i±iin, rectangular and skew h^rid laminated 
composite plates with and without mixed boundary 
conditions have been obtained. They found that hybrid 
laminates provide stiffnesses intermediate to that of 
single material laminates and that the first natural 
frequency of the hybrid laminate can be made higher than 
the stiff er single material laminate. 

1.2.3 Second Order Shear Deformation Theory 

29 

Whitney and Sun proposed a second order shear 
deformation theory ^^ere in two terms were considered in 
the series form of the displacement function as follows. 
The effect of thickness on^displacement was also 
considered. 

u = «Q(x,y,t) + {x,y,t) + ^ 0 ^ (x,y,t> 

V = + Zf^ (K»y»t) + 0y (x,y,t) 

w = % (x,y,ty + Zf (x,y,t) 

They obtained the fundamental frequency for various 
laminate layups with simply supported conditions. All 
these results were compared with those obtained from 
classical anisotropic theory and the exact results of 
elasticity solution. It was found that ihis higher order 
theory yields results comparable with those of elasticity 
solution. 
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1 •2*4 Higher Order Shear Deformation Theory 
26 

Lo et al developed a higher order theory for 

analysing the plate deformations* They in addition 

to the higher order terms considered by Whitney and 
2i6 

Sun , considered a third order term in the displacement 
function for inplane displacements and a second order ' 

term for the lateral displacement function ’w’* The 
form of this displacement field is as follows: 

u = UQ(x,y,t) + Z. <^(x,y,t) + jj(x,y, t> 

+ 0^ (x,y,t) 

V = %(x,y,t) + Z.fy(x,y,t) + ^?y(x,y,t) + ^0y(x,y,t) j 

I 

w = w (x,y,t) + Z.f (x,y,t) + (x,y,t> j 

o Z - Z I 

They used this higher order theory for the static analysis ! 

I 

of laminated composite plates and obtained results 

comparable with the exact results of elasticity solution* | 

27,28 i 

Similarly Karaal and Durvasula presented a [ 

i 

higher order shear deformation theory for -Uie dynamic analysis 

, 

of rectangular and skew composite plates with simply supported j 

I 

and clamped boundary conditions* Here they considered a 
displacement field such that the condition of zero transverse 
shear stress at the plate surfaces was satisfied* The 
displacement field considered is as follows. They did not 
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consider the effect of thickness on lateral displacement* 
u = Uo(x,y,t) + fl(Z).fj^+ f2(Z) 

V = + flCz).!' y + f2(Z) 1^ 

w = vr(x,y^t) 
where 

f 1 (z) = j • Z — -- ^ TT • 2 ^ 

^ 3h"^ 

f2(Z) = i . Z - ^ . Z^ 

^ 3h^ 

Rayleigh-Ritz techni®i» is employed for obtaining the 

of 

natural frequencies and the mode shapes_/simply supported 
and clamped, rectangular and skew composite plates. 

Results so obtained are compared with classical lamination 
theory and other higher order theories in which Finite 
Element technique was used. It was found that the 
results were comparable and that ilie effect of thickness 
on lateral displacement was negligible* 

29 

Reddy also proposed a higher order theory which 
is unique in comparison to the earlier meaiaoned theories* 

Even though, iliis displacement field very much resenfcles 

Q.W ,9S 

■fee one proposed by Kamal and Durvasula ^ , it differs 

in one respect* All these earlier mentioned theories 
were based on assiraed displacement field and there was no 
theoretical backup for feese displacement fields* Although 
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Kamal and Durvasula^^*^® satisfy the condition of zero 

transverse shear at the laminate surfaces, they assumed 

a displacement function to satisfy this condition. But 
29 

Reddy, using the condition of zero transverse shear at 
laminate surfaces derived the displacement field. The 
displacement field is as below. A detailed explanation 
of this model is discussed in adapter 2. 

u = UQ(x,y,t) + fl(z). p^ 4 f2(Z) . H 

V = Vo(x,y,t) 4 f1 (Z) . 4 f2(Z) o 

w = Wq (x,y,t) 

where 

■ \ 

fl(z) = z - 

3h 

f2(Z) = - . Z^ 

3h 

It is observed here that this resembles the displacement 

field assumed in Ref. 27 and 28. He Justified the 

assumption that 'w* is not a function of 'Z', since the 

is 2 

transverse normal stress /)f the order of (h/a) times 
the inplane normal stress. Navier approach has been 
employed to obtain exact solutions of -ttie equations 
obtained by using the principle of virtual displacement for 
symmetric cross-ply simply supported laminates. Results 
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for static analysis have been obtained and compared viith 
three-dimensional elasticity solutions of Pagano with 
simply supported edge conditions under different loading 
conditions. 

1 .3 Objective and Scope of the Present Work 

The foregoing literature review clearly Indicates 
that very few investigations have been carried out on 
the static and dynamic analysis of hybrid laminates and 
most of the studies either neglected the transverse shear 
deformation and the rotatory inertia effects, or a first 
order shear deformation theory was used. In this thesis 
an attempt has been made to analyse the free vibration 
response of hybrid laminated composite plates. Typical 
results have been presented for simply supported 
boundary conditions. The higher order theory proposed by 
Reddy^' has been employed to take into account transverse 
shear and rotatory inertia effects. Rayleigh-Ri 
technique has been used to <±>tain the natural frequencies 
and mode shapes. The results obtained have been compared 
with available results in iiie literature. To check 
the accuracy of the technique, study has been made for 
isotropic plates and regular composites* In all these 
cases the effect of length to thickness ratio, aspect 
ratio, nuBPber of layers in ihe laminate for symmetric 
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and antisymmetric angle-plies with dif f erent fibre 
orientations and for symmetric cross-ply laminates 
have been obtainedo The effect of position of a 
particular lamina in a laminate also has been 
analyieed. As a special case results for a two layer 
isotropic hybrid laminate has been obtained. This has 
been accomplished by a programme developed on DEC-1090 
System using standard ligen-value routine and reported in 
subsequent chapters. 

In Oiapter 2 mathematical formulation has been 
discussed in detail. Simply supported boundary condition 
and other simplifications are also presented. 

In Chapter 3, results have been discussed and at 
the end scope for futurework with conclusion has been 
included. 
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CHAPTER 2 

FORMJLATig] AND ANALYSIS 


2»1 Shear Deformation Theory 
2.1 .1 Introduction 

The classical lamination theory is based on the 
Kirchoff s hypothesis which states that, 

o Normals to the raidplane before deformation remain 

straight and normal after deformationo 

• Stresses acting in the XY-plane dominate the plate 
behaviour, i»e», a x and are negligible. 

o Displacements u,v and w are small in comparison to 

the plate thickness. 

« Strains e„, e„ andy are small when compared 

X y xy 

to unity. 

• Rotatory inertia effects are neglected. 

It has been well established from experiments that the 
above assumptions result in the under prediction of 
deflection aid over prediction of natural f reque rwi es . 

This is due to the fact that the traasverse shear strain 
and rotatory inertia effeetts are neglected in the classical 
lamination theory. The transverse shear stresses although ‘ 
small cannot be neglected and can be calculated by 



17 


integrating equilibrium equations. A nondimensional 
analysis has shown that the transverse shear stresses 
are of the order of (h/a) times the inplane stresses 
while the transverse normal stresses are of the order 

\ 

of (h/a) times the inplane stresses, where 'h' is 
the thickness of the plate, 'a* is the length of the 
plate. From this we can assume, 

(i) The transverse shear stresses can not be neglected 
for all sizes of plates. 

(ii) Transverse normal stresses can be neglected since 

the order of magnitude of these stresses in comparison 
to inplane stresses is very small*. 

The errors in buckling load, natural frequency etc. are 
quite high for plates made of composite materials whose 
ratio of elastic modulus to shear modulus is very 
large (e.g. Graphite-epoxy, Boron-epoxy). Hence, for 
analysing the structural behaviour of composite plates 
one has to consider a theory which takes into account 
the transverse shear effects. 

2.1.2 Higher Order Shear Deformation Theory 

The higher order liieories have been developed 
using the kinematic displacement field defined by 
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u(x,y,z,t) = u^(x,y,t) + 2 p^(x,y,t) + §jj(x,y,t) 

O 

4* z 

v(x,y,z,t) = VQ(x,y,t) + 2 py(x,y,t) + z^§y(x,y, t) 

+ 2^ 5 y(x,y, t)+* •• • 

w(x,y,z,t) = WQ(x,y,t) + zp 2(x,y,t) + 

+ 2^?2 (x,y,t) + .« .0 . 

(2.1) 

It can be seen that with each power of the thickness 
coordinate an additional unknown is introduced into each 
component of total displacement. The equations derived 
from the Hamilton’s principle would become difficult to 
handle and for many unknowns p^iysical interpretation can- 
not easily be given. Also these theories do not satisfy 
the boundary condition that the transverse shear stresses 
should be zero at plate surfaces. For these reasons, in 
this work a refined higher order theory proposed by 
Reddy has been considered. This theory is also based on 
the displacement field given by equation (2.I) but the 
specific form is chosen to satisfy the boundary conditions 

''xy = i |) = ''yz =± |) = 0 

This choice eliminates the need for any shear correction 
factor and third order distortion of liie normal to the 
midplane is allowed. 
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We consider a displacement field wherein the dis- 
placements along x- and y-directions are ®xpaade4 a*,, 
cubic polinomials of the thickness coordinate and 
displacement in the tramsverse direction is assumed to 
be constant through thickness due to -Uie reasons mentioned 
in Section 2.1 *1 

2 , 

u(x,y,z,t) = UQ(x,y,t) + z p j^(x,y, t) + z g^(x,y,t) 

+ z^ ? 3j(x,y,t) 

v(x,y,z,t) = V (x,y,t) + zp (x,y,t) + z^4 (x,y,t) 

0 I y 

w(x,y,z,t) = WQ(x,y,t) (2*2) 


where u^, and w« are the raidplane displacements, 
o o o 

S and B are the rotations of the normal to the midplane 
x y 

and the functions k and have to be 

determined from the condition the shear stresses x 

and X be zero on the top andJaettom surfaces of the plate#, 
yz 

For orthotropic plates this condition is equivalent to the 
requirement that the corresponding strains be zero on 
these surfaces* We have. 



+ 
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yz 


9 V 

8z 


3 w 

ay 


= ,, + 22 £, + 3z^? 


3Wc 

9 z 


(2.3) 


Using the condition that these strains be zero at 
z = + h/2, we get 


= P + h S + 2^ 

X X 4 

fx - 

3 w 

Q 

3 X 

at z = + h/2 

= p . h S + ^ ? 

X X 4 X 

& X 

at 

z = -h/2 


From these two relations, we get 


3w. 


0 and ^ + ^ ] 


(2.4a) 


Similarly using the relation for Y , we get 

y ^ 


S = 0 and ~ ^ ^ 


3h 


3^0 

ay 


+ T-^ 3 


(2.4b) 


Substituting relations (2.4a) and (2.4b) in equation 
(2. 2), the displacement field becomes 

3w 

« = P X * jf • 

V = Vq ■^ Py • f^(2) + * f2(z)l 
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w = w 

0 


where 

fl(z) 

= z - 0 z^ 

3h 

f2(z) 

li 

1 

• 

f 1 (z) 

il 

0 

• 

N 

1 

f2(z) 

= -04. Z? 

where Cl 

= 1 and C2 = 



(2o5) 


(2o6) 


The displacement field used by Kamal and 
97 98 

Durvasula * shows that only the form of fl(z) and 
f2(z) changes since they assumed this form to satisfy 
the condition that t and t to be zero at the surfaces 
while Reddy obtained the form using these conditions 
as shown above. f1(z) and f2(z) found in Ref *27 and 
2% are as follows S 


% 5 _ 5 _3 

xUz) = j z - — w z 

^ 3h 

f2(z) = T ^ 

^ 3h 

Thus it can be seen that even though cubic variation 
of the inplcffie displacements are considered, the number of 
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dependent variables remains same as in first order 
shear deformation theory. 


2.2 Strain-Displacement Relations 

The strains corresponding to displacement field 
defined by eq. (2.5) are 


^ X 5x 

— JL 

~ SIT 

fl(z) • 

IS 

ax 

f2(z) . 

3 

a x^ 

- 

^ ay 

= 

ay 

fl(z) • 

ay 

to 

N 

• 

-2 

ay 


xy 


a u . 3v 

sy ax 


a u 


3V, 


Y = X + 

^ xy ay 


+ fl(z) .[ 


iEx + + 2.f2(z). 


ay 


ax 


xz 


= Jji 
a z 


a w 

o 

ax 


JW 

a xay 


] 


= fl '(z) . P „ +, f2'(z) . 


I ^ 

' o 

a X 


a w 


a X 


aw 


yz 


fi‘(z).Py 4 


^ [1 4 f2’(z) ] 

" [ 1 4 f2'(z)l 


(2.7) 
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Substituting for fl(z) and f2(z) from eq. (2*6), 


we 

get 

the strains as 

€ 

X 

= 1 

hi ^ ^ 


= 

21 + * [•'21 

Sz 

== 

• ‘=42 

Y 

xz 

= t 

•51 ^ ^^•‘'52 

^ xy 

= 

% * ^C'ei 


= 

3 u 

it ’ "ii =«’ 

®21 


3 V 

■77”=2i =='• 

^41 

= 

a V 
Ty 


= 

d w 

(Cl . Pjj Tt ■ 


: k.. = -(C2. + C. 


3 V 


. 2 

3 w V 

“T ) 


3u 3v^ 

q . O » ir 

1y 35r * ' 


36 

+ 


= - I G2 ( 


^ V 2 C4 ^ 1 

77 + it ' + 2.04 . ] 


( 2 . 8 ) 
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2o3 Stress-strain Relations 

A laminated plate consists of a number of orthotropic 
plies -with perfect bonding. Each ply will have principal 
orthotropic axes at an arbitrary orientation to the 
reference axes of the plate. If the state of plane stress 
is assumed, then inplane stress-strain relations can be 
written as 

k 

■<ri. 

, fT 

Xlt) 


°11 

CM 

a 

1 

0 


>L ] 

°12 

^2 

0 

i 


1 

0 

0 

%6 


I'ul 


and the transverse shear stress-strain relations can 
writtom as 



be 


(2.9) 


where, 

^11 “ 

^22 “ “^LT^TL^ 

^12 ^ ^LT’ LT*TL^ 


I 



25 


Q 


"66 

'■44 

^55 

k 


= G, 


IT 
^TZ 

\z 

kth layer, 


( 2 . 10 ) 


These coefficients are known as reduced stiffnesses. E,G, 
y are the Young's modulus, shear modulus and Poisson's 
ratio respectively. The subscripts L and T are parallel 
and perpendicular to the principal orthotropic axes of 
each ply. 

The stress-strain relations for laminated plate 
consisting of 'n ' nuniber of layers at arbitrary orientation 
to the reference axes are obtained using transformed 
stiffness matrix. 


X 


Ql1 

S2 

^16 

0 

0 

/ 

4 ^ 

^X 

^ y 


S2 

^22 

^26 

0 

0 

I 

fiy 

V ^ 5 Cy 


^16 

^26 

^66 

0 

0 

' 

f I 

xy 

yz 


0 

0 

0 

Q44 

^5 


f 

yz 

, *’^xz 

^ * 


0 

0 

0 

^45 

%5 


f 

’xz 

V J 


( 2 . 11 ) 

where «re given by 
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Ql1 = +2(Q^2 + ^22^"^ 


‘22 = Si + ^<S2 ■" 2^66 


S2 “ ^Sl S2 ” S2^^'^ 


*^16 ^Sl ” S2 “ ■'■ ^S2 “ S2 '*' 


^26 ~ ^S 1 “ S2 “ ^^12 ” S2 


^66 “ ^Sl S 2 " ^S2 “ + 5 “^) 


Q44 — ^44*^ + Ss*^ 


— P44 + Q ‘C^ 


55 


Q. 


'45 


^S4 - 


(2aJ2) 


where, C = Cos(©) and S = Sin(9) 


and ’©' is the angle between the reference axis and 
principal or thotropy axis as shown in Fignre 2.1* 


2 *4 Plate Constitutive Relations 

The stresses in a layered composite vary from 
layer to layer* Hence it calls for a simpler bwt statically 
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equivalent system of force and moment system* acting on 

the cross section of the laminate vhich can be used for 

Such a system 

exact analysis of these structures^as been shown in 
Fig. 2 . 2 (a) and Fig. 2.2(b). The resultant forces and 
moments acting on the laminate cross section are defined 
as follows! 



"v 

N ) 
XY 

h/2 

= / 

-h/2 

(a , 0 ,T )dz 

X* Y 


(Qx 

.Qy) 

h/2 

= / 

-h/2 

f ^ ) dz 

xz* Y^ 


M 

Y* 

“xy> 

h/2 
= / 

-h/2 

( 0 ^ ) .z dz 

X* Y ^Y 


Py» 

fxy) 

h/2 

= / 

-h/2 

(cr ,0 ,T ). z^ dz 
x’ y • xy 

(R^. 



h/2 

= / 

-h/2 

(t ,T ). Z^ dz 
xz’ yz 


(2.13) 


Using the relations, Eqs . ( 2 .II) and (2.13), we get the 
modified constitutive relations which can be used in 
the buckling^vibration of plates and shells. 


Uy] 




^44 


■°44 °45'f 

1 

' 41 ] 

1 

l°x] 




J^45 ^5 


L°45 “ 55 ] 


® 51 ) 

/ > 
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^ '' y i 




^4t °45l 
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F F 

^44 ^45 

1 

■i 

42 
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> °55_ 


^45 ^55 
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\ 

^2^ 
V J 

m 


(2.144 
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N. 


N 


xy 








N 2 ^16 


^2 ^6 
^6 ^66 


B 


11 


B 


B 


12 

16 


®12 ®16 
®22 ®26 


®26 ®66 


D 


11 


D 


12 


D 


16 


^12 ^16 
^22 ^26 
^26 ^66 


1 1 


12 


16 


®12 

— 

®16 

IP— 

^11 

"12 

=16] 

fre \] 
11 

®22 

®26 

®12 

^22 

>^6 

^21 

®26 

®66_ 

^16 

u 

®25 

=6^ 

1^61^ 

“ l 2 

“id" 

^11 

52 

=16] 

(k ^ 
^11 

^22 

°26 

^2 

==22 

=26 

hv 

°26 

°66. 

Jl 6 

==26 

==66 

ki 

^12 

^16' 

»n 

52 

=16] 

(52) 

’^22 

==26 

^12 

=22 

=26 

521 

^26 

==66. 

” l 6 

=26 

=66] 

( 52 ) 


’where 


N k _ 


2 3 4 6v^ 

yZ yZ j ciz 


for i = 1 ,2,6 


N 


(Ajj.Dij, [Q.^] (I,z2, z‘’)dz 

for i,j = 4,5 {2.15) 


2.5 Strain Energy and Kinetic Energy Relations 

The strain energy of a rectangular plate can be 
expressed in terms of the stresses and strains as 
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U = i 


a b 

/ / 

o o 


Ul 

z 

k=1 


/ 

^k-1 


^ 4*' T 

X ^ y y xy xy xz xz 
+ 'fyZ Yy2> '^y 


(2.16) 

Using the equations (2.11) and (2.14) we get the strain 
energy relation for the rectangular laminated plate as below! 


a b 


3 u 
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3u 
( 

3x 


2 - { 

3v 


3v 


o 3 u 

+ 2A ^ 

ay ^'^IS 3x 


^ o 3v au a V 

("T^ 7T) + (—^)^ + 2A_ — (— ^ 4 — ^ ) 

'ay ax' ^2 ZY ' ^ a y a y ax' 




+ 2 A 45 (c^p ^ c, |H ^ c, * |^.|^) 
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a u 
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X 10 
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a V 

^ V 


9P _ 

au 8 v 

^ ^®26 <-, 

0 

ay 

(-i_-S 4 

' ay 


-sf V ^ 


a 

ap^ 

36 

a V 

dp ap 

®66 

0 

a y 


•" 2ci gri 

<■77 “5^ 

9 

2 

2 

ap 

ap 


+ c, D,, ( 

X\^ 

“57^ 

4 2c^ 

^12 ax 

-_Jt 

ay 


0 

ap 

ap 

ap 

0 

*P 0 

" 2c; 0,6 

X 

ax 

(«_x 
a y 

Tx> - 

°22 '1 




ap 

ap 



2 

4 2c^ ^26 

Zi 
a y 

<T^ ^ 

-^) 
a x' 




-6D^4 3 y ^ y By 

/3w\2x 
^4 ^ay^ ^ 


-6 D 45 ( 2 c^ ^2^x^y‘^^2Pxay '^^l‘^4^x ay 


av 


, ^ _ w rr ^ O T n » O ^ ^ ^ w V \ 

^ ^2 ^ y 51c ^1 ^4 ^ y a ■“ 


3 J£ 

X 


a V a v ^ 
'4 a X a y^ 


%5 ^2^ X ^2 ^ X ax ■‘' ^1 ^4 ^ X ale ^4 ^Ix^ ^ 

2p ' ^ 3 P ^P 


- 2E41 (C2 


ay 

au 


O ^X ^ ^ 3 % a ^ 


31 


a u 

- 2Ei 2 (=2-1-1 • Ty 


a a^w 

* * ay2 

a V ap 

*2 ay -an * 

c 

4 ay 

a^v^ 

a u ¥ 

-2Ei 6 (<=2 TJ ^"jy ■' 

a x' 


a u a ^w 

^ ^^4 a X a 36 y 2 

aa au a v 

Z2L (-Ji + -J 
ax ay a 

a\ a 

'4 j-T- (T? " 

av 

-t5 » 


.2= , (c . 

^-22 '=2 a y ay 

av 

c — 

4 ay 

—2-) 

ay 


- 2^6 <'2 


ap. 

Vx 


a V 

4 ' iy a xay 


V 



ap 

d U 

av 


a 2^, 

au 

3 

4 C 2 

— ( 
a y 

-r^ 4 

ay 

0 ^ 

^=4 

a w / 
”T ^ 
a y 

T7 ^ 

0 

a X 

u 

-66 

( 2 C 2 

a u 

i-a 4 

ay ay 

ax^ 

4 40^ 

a u 

0 

a y 

2 

a w 

a 36 y 

-* 2C2 

av 

a 

ax 

a^ 

<-a7 ^ 

dp 

ax'' 

4 4c 

9vo 
4 a>^ 

2 

a w 

a xay 

) 

1 

(CiC^ 

(!Li)^ 

> V ax'' 

4 c 

1^=4 

d| 

X 

: 

2 

a V 

y , 1 ■ 

3 



-2Fi2 (2c, C 2 -3^ -sf |i<=4 j 2 

y ^ ^ \ 

+ 2C4C, ' — o' * 

^ ^ ^ ax 



32 


-2^6<2'=i'=2 ^ » « 


2 

^ V 

3 36 y 


~2 JS 36 

- ' 1'4 4 -^)) 


-2F22 (c,C2 (- t ^)2 + c,c^ -5^ ^) 


a? a# 


- 2 F 26 (2c, C 2 - 3 ^ + 2 c, 


ay a xay 




^ c,c, 3 (-^ . 


91^44 ('iPy *'4 (|7>^ 4 2'=2'4P 


'4 '3y^ 


c ft 

2^4^ y 3 y^ 


I8F45 (c2 -P P y + p 


2^4 P X 3 y ''^ ^ 2^4 ^ y 3 X 


.2 ^ av^ 

^4 ax 3y^ 


A 2 rt 2 ^ f o\f\^ . ^ ^ ft c w' 

"^55 4 2c2^4 Px 3 x^ 


.2 /3w\2 


ap,,. ae^2 

'2?66 (‘=i'=2 ) 


ap ap 2 ?p ^ 

42c, C^ 4^ -df + 2C,C^ (-J-^ 4 -^)) 


'1^2 ay ax 


4H„ (c^ ^2C2C, - 


, 2 

X , 2 ^ 

O X 



33 
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(2.17) 


Similarly the kinetic energy of the plate can be written 


as 


T = ST 


a b h/2 
Iff P 
o o -h/2 


(k) 
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(k) 

where P ' ^ is the mass density of the k layer. Now 
using the relations for the displacements u, v and w 
from equation (2*5) we get the kinetic energy relation 
for the plate as followst 



where P, q, and I are the normal, normal-rotatory and | 

rotatory inertia coefficients given as 

P, al, 02, 11. 12. 13 = S [1 ,fl(z),f 2 ( 2 ),f 1 ^{ 2 ) 

k=1 

f2^(2),fl(z).f2(22 dz 

(-. 20 ) 
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2o6 Rayleiqh-Ritz Method 

Rayleigh's energy method is used to obtain the 
fundamental frequency parameter of a system for which 
solution of the eigen-value problem is difficult to 
obtain. This method gives a higher estimate of the 
fundamental frequency parameter. The closeness of 
this value depends on the degree of resentolance of 
the trial function to the natural mode. 

The trial functions which satisfy all the 
boundary conditions of the problem and are differentiable 
as many times as the degree of the system are known as 
'comparison functions', while those functions that satisfy 
all the geometric boundary conditions and is differentiable 
half as many times as the degree of the system are known 
as 'admissible functions'. Hence the trial functions 
selected in Rayleigh-Ritz method must belong at least to 
the class of 'admissible functions'. The trial functions 
chosen in this method will be of the following forms 

M N 

u(x,y) = S 2 a. . C(x). f .(y) (2.21) 

i=1 j=:l ^ J 

where • are the unknown coefficients to be determined. 

0. (x) and * .-(y) are admissible functims. These functions 

mim ^ 

0£(x) and ^jCy) knowh as generating sets and the 
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coefficients are determined so that u(x,y) resembles 
as close as possible to the natural modes. This is 
mathematically equivalent to seeking those values of 
a^j for which Rayleigh’s quotient is stationary© 

The Rayleigh's quotient can be expressed as 


R(u) 



N(a^^),l.J=1.2,...n 
D( 1 f 3 — 1 * • o .n 


vjhere. 


V__y is the maximum strain energy 

Xu 

* 

T is the reference kinetic ®nergy 

^ is the fundamental frequency. 


It can be seen that the admissible functions can be integrated 
over the spatial domain involved in strain energy and kinetic 
energy thus rendering Rayleigh’s quotient as a function of 
the coefficients a^j (i, j = 1,2,...n). The quotient will 
have stationary value if the variation vanishes 


3 R 



3 

<3L ^ 



n{- 


3d 


30 / 


rs 


= 0, for r,s=1,2,«.on 


2 * 

Denoting the value of » associated with the stationary 
value of the quotient by A , we get 


3 N 


3 D 


rs 


= 0 


(2 . 22 ) 


Here 'N' denotes ilie strain energy and 'D* demtes the 
kinetic energy© 
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Equation (2.22) is used to obtain a general solution 
to the vibration problem of a laminated plate. We assume 
a double series form of admissible function for the displace- 
ments iIq* ^o* ^o rotations and as follows: 


u. 
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= 2 S 
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mn ym yn 


(2.23) 


where i and h are the normalised values of x and y given 
by g = x/a and n = y/b. 0 ^ 3 C£ ) , 

of ( t) 0 ( S) are assumed as trigonometric functi^s for 

* yf 

the inplane displacments (u^ and v^), lateral displacement 
(w ) and rotations ^ jj and P^) respectively corresponding 
to the m mode satisfying the boundary conditions at 
I =0 and g = 1* Similarly the ? -functions corresponds 

to the n?^ mode satisfying the boundary conditions at 
T} = 0 and =1* Thus depending on the boundary condition 
in each of the g i} directions one can choose the 

, -.1 
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corresponding admissible functions that will satisfy 
the boundary conditions* The strain energy and kinetic 
energy expressions will take the following form after 
equations (2*23) are substituted in relations (2.17) 
and (2*19). 
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2 


+ 2.B11 . {C1. a . ^ ^ 


rs xr 


0' (£) 


xs 


+ 2.B12 .tC1 . a^ . * ^rs’^yr^^ ^ 



+ cUa^A Gjjjn- 0 ^ib(§) ‘ ^rs* ^ 


4* 


2.BI6. tC1 •[ a^A. C . -^yni^ £^*’®un^ ^ " ^rs’^xr^®^ *'*xs^’^^ 


mn 


+ C^„. (£) • ^ra- •»ys<’’’ 

- A- >«mn- 


4- a 
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+ a' 


?. K.. (5). *,„{-!) . 0^3- 0^3- 


^xm 


+ 2B22i C1 .a^.??o Gjjj^.0^jn( l) • ys^ ^ ^ 


.826. { C1 { a^. ^ 2^vin*^vm^^ ^ * ^rs*^xr^ '*xs^ 


. X. G_.0„JI) 


mn vm 


+ a •\ 


.X^ . Y... 0.™(£) . *i.(T.)- 0^3- -’-is'”’ 


mn ym 


+ a®. X. Y„„. 0^(S) °rs- 


+ 2.B66. t C1 ^ • C n-^untf) • ' ^rs'^’xr^ ® 


.'*'’ (n) 
xs' “ 


- a^X. C3„- • ''rs- 


^ 2 . p (a* (S') « (ti) . X . 0^Jt) 

+ a .X. G^J^* "^un' ' rs xr xs 


. a^G,„. 0J„( a • 

^ D,1.{ C1^ a^ . X^. 0ij£) ».„("> • ’^s- •*xs‘’» » 

+ 2.D12t C1" . a2.x. X,„. 0i„U) * ^ra-V‘'^’ ya'”^ ’ 


+ 2-D16 .t C1^ .[a^.X* Xgy^* 


. iSji) =^rs- 


+ a^. X„„. 0y« •»xn‘’> • '^ra* * ’ya^’'> 1’ 
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+ D22 .{ Cl2 . . y^„. 0y„,(£ ) • »^(n) . 1)1 

+ 2 .m 6 .t Ci2 { Y„„. 0 ^( s) .j^n) . ).*'^(n) 

■" \n' '^rs* 

- 6.D44.C C1.02.a'*. 0^(i ) «^(’)) . Y_.. 0._ ( S) .»„( 1) 


rs yn 


ys 


+ C2.a .X. Y„„. 0^(1) .*^(11) . 0^,.A_(6) .?,'„(„) 


yn 


rs wr 


ws 


+ C4.[ Cl.a .X. .W^(n) . O_.0_(§) 


yn. 


rs wr' 


ws 


2 2 


^ • °m- a .?i,(Tl) . 0^^. 0^(S) . *.Jn )]1 


rs '^•wT' 


ws 


- 6.D45 .C 2C1.C2.a 0,^(£ ) 


rs rs 


yS 


+ ^ * °rs*^wa^^^ CUC4) 

+ 2.a^?u o^„. 0 ^{%) .*^(d) . o^^.0^(s) ar;^(ii) . C4l 

- 6 •D 55 .^ 01 .C 2 .a^ X^„. 0 ^„(£) .*^„C 7 ,) . X_ 0 ( 6 ) (n) 


+ C4. 

2 

a • 

o«« 0’ (i) f ^ 

(nj .0^ .0* 


MU Win wn 

rs wr 

. 3 

+ a • 

\n- 

0 (t ) • (^) 

. O^^.0» (1) 
rs 1 ® 


WS 


ws 
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+ D 66 o{Cl' 


r 2 2 

^ i 


X 


mn 




X„.^„C£) 


•lyn) 


"■’ ^mn- •®Yn^’’> • « . V 'l) ^ 

+ 2. Cl^ . a^. .5 y „) , P’yr'^S) ^'^(<1) 

- 2.E11 .t C2. a2. C_^. £) ^)y f) 5 ^ 3 ( 11 ) 

+ C 4 . a. c^„. 0 '^(|) . 0^^. ^ys ) .,^^(n) 1 


-2.E12 .{C2. a^. X. C^. .I^nCri ). la^^Cf) ^^>1) 

+ C4. a.X^ C^. )»^^(£) .I^^C’I) . 0^^. 0^^(£) .IJ^^Cn) 

+ 02 . a\. G^. 0 ^C£) . sy „) . 0 y£) .T^^Cn) 

■*■ '=''• °mn' •* • Ors* > • » ws^” > ’ 

- 2 .E 16 .tC 2 { a^ „ X, Cji^. 0 >;^(£) .1 . Xj,^. ) 8 jjj ,(£) 

+ '^mn- P'u=^5> -funC) • \s- • ® ys^’' > 

+ a^.X . X^. 0y £) .T^(ri) . p^^£ ) 

>W- • °ra- 

**' 0-^,2,aX.. Cjjjjj. £) •fynC ll) • Oj.3* ?^g(^) 
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°mn- • ^rs- > •» is^l ) 

+ ^ • °mn- • S„. /»'„(£) 

2E22 , t 02. . x2. G_^^. 0J, {) . w^Ct)) . £) * ) 

+ C4 . a.x3 . ,i;^(„) J 

-2.E26 , { 02. [a^A ^ . G^„. £) . ( >1) . 0„(£);j'i^( 

+ '^- • '^rs- ^'yj- 0 •*ys<*’> 

^ ^ • ^rs* ^ 

+ a2.x. Y„. 0^(£) . ] 

+ 04. [ 2a. . G^. £) ., ^Cn ) .0^3. p-^E) jjy ,) 

v| 

+ a*X * Ojjin* • ^v«n^ * ^rs* ^ 

+ °mn* •’*' ^ « °rs* 

— 2,E66,i C2*[ a ,X ^ Am * A^ujCS ) • ^rs* 

» 

+ a .•x* ^lan* ^xxJ' un^ ^ • \s* ^ • 'f ygC'*!) 
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• °mn" •* y” ^ 

■" • vn^”’ - ■ • "^rs- ^ 

+ 2. C4.[ a.x2. C^. JS^CJ) . ^ • •^sC'l > 

* ‘\m- • *y ’l) . Ors* ^'wr^ •* ^ ” > 3* 

-2.F11 .tC,.[C2. a2. p^(j) .x,3.^y£) 

+ V- py f ) • lyi) . o^B- »* ws^” > ]» 

-2.F12 . { Cl . [ 2.C2.a2.V. X^. J) .,^(„ ) . 

+ C4. a;, 2. J-yi) . 0^^. £) .i y ) 

+ C4. aA. la^(J) . ^(„) . 0 ^^. J) .J„g( Ti) ]1 

- 2.F16 . (Cl. [2C2.t a2.K. X^. t) . * ^Cn ) . £) .,y 

+ \n-<®y «) - ’■y’l) • ^ra- ^yr« > ’ 

+ 2C4, aXi X^. £)i„(£) .f^(n) . 0,.^. (5;^£) .Ti^Cn) 

+ C4 . (a. X. q,„. ^;^C£) .ijn) . 0^^(£) ;t y r,) 

+ "• °mn- • y ”) • \a- •»ys^’'> ’1 
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-2.F22.tC1. [C2. a^. J T] 

'^mn- • °ra- l^wr^ ^ 

-2.F26 .t C1 ^ 2.C2.taV. Y^. i) y ^(„) . X^^(S) .,.^(„) 

•" ^mn- • ^rs- ^;r^£ ) . I yjCtl ) 1 

+ C 4 t 2 a. . Y^. IS^CI) .y y „) . o^^. (£) ^^3(11) 

+ a.X^ . O,^. PJi ) .1'^t.r>) . X^3. P^J.S) . S) 

‘^n- • Yj-g. P'yjXS) ^ysC’l) n 

+ 9. F44.[ 02^. a''. Y^. p^it) . Y^3. P^J. i) ’n) 

2 I 

+ c^ . a ^ . 0^„ Pj.i) .yj^(a) , 0^3. P^j^t) a) 

H- 2.C2.C4. a^.X . Y3,„. P^ii ) . 0^3. VC « * Vs ^ ’ 3 

+ 18 F45„[ C^ . a^ . V., ) J'xr.C*') • l^rs- V<* ^ » taC" > 

4 C2. C4.t a^^.. V «>„(£) . Ors- ^wiC e •» is^l ) 

■*■ a • ^jnn* ^yifiCS ) y yjjCi>1 ) • ®rs* ^*wrC£3 •* y,gC ’1)1 
+ li. a X » Oj^n. Vi(£) J ^C’)) . 0 ^, 3 . 03^( i) ,l ) ] 
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+ 9.F55,. [C|. .J^^CTi 

+ ^4 * ® • ®mn* ^ ^ ^ * ^rs* 

+ 2.C2.C4. a3. X^.PJi) .^^(r,) . 0^^. P'^CS) .1 ^3(11)] 

- 2,F66.. !C1.C2.[ a2.x2^ x„,„. n„(«) • >^rs -^xr^ > 

+ '^mn- • vn^l*) • Vr3. |Zi'^(S) .*^3(111 

+ 2. Cl .02. a^.x . x„^. 0^(1) ^^(r, ) 

+ 2, Cl .04.[ ax® . 03,„. 0 ;^ (S) S)* yn ) 

+ a X. C^n. 0'„(l) .Si„(il) . Y^ 3 . 0^j.{S) , JygCn)]! 

+ H11 .t C® a2. x„„. 0^ (E) .I^Cn) . 0y£ 

+ =4- °mn- pys ) • ’wn^”) • Opa- KJ- •’ »=(’> ) 

+ 2.,C2.C4. a,. Xj 3 „. 0yE) .S^Cl) . 0 ^ 3 . ^-^^CE ) ,.l'„ 3 C’l) J 
+ 2. H12 .tc| ,. a^X. X3,„. 0'j.i) .w ^(r,) ,. Y^^. 

+ G2.C4, [ a> • .f^nCn) * O^^g, ) *fj,g(tl) 

* ®rs* ws^ '^ ) 3 
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°4 • 0„n- . 0^3- 

+ 2. H16, CC^ . [a2 , x. p-^CS) . ) . X^^.P^^CS) 

.fis( n) 

V>* • '^rs- P'yr‘^) •» ys^” > 3 

+ C2.C4[ 2.aX. X„„. p-^(| ) . T^„(.,) . 0^^. fi ’ )., 

* ^ V- .f„(t,) . X^3. P^J.^) .t y r, ) 

■*■ °mn* ^ • ^a* • *ys^ ] 

+ 20^ .X . 0^„. 0"^( t) .S^(T)) . 0^3. p'^^CI) . fy T, ) I 

+ H22 . tC® . a2.x®. Y^„, Q^. P^(£) .»;„(t)). Y^3.Py^(£).»yT 

+ C4 , X • ®j]jn* ^vjm^ ^ • ^rs* *'®' ws^ 

+ 2.C2.C4. a.X® . Y^n. 0^(£) f ^ H ) , 0^3. P^(S ) ar;„(ll) J 

+ 2.H26. [c| ,[ a®.x Y^^ti. |5ynC£) . ^rs -’is^ 

+ 4 • X* ^uin* ^ym^ yn^*^ ^ • ^rs“ ^ 

•I* C2 2,^ ,'B^X. '' • • ®rs* ^ 



<N CM 
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+ a. X . Q 


mn' 


• °nin- 


+ 2. c| . 


+ H66 f 4, C4 . 



® ^ * ^rs* 

'“«(«) •» P;,r(£) . ] 

°mn- *;:„(il ) . 0^^. la^j) . i; 3 Ct,) 

• °mn- ^’wm •5’v*^’' > • O^s* 

• \n* ® ^rs* *'*'xs^''^^ 


+ '^mti- • V”’ '^rs- ^'yr(» •»Yn(’l) 3 

+ 2, C 2 • a ,. x, ^nin* ^xm^ * '^rs* ^ 

+ 4. C2.04. ax^. |8yi) «;„ (n) . X^3..0^^({).f;^^(D) 

+ 4.C2.C4. ax . 0^„. 0;^ (|) ,;^C„ ). P'^p) Sy^(tlJ]d£d 



(2.24) 
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The kinetic energy given by equation (2*19) now 
takes the follov/ing forms 

(P, q1,q2, 11,12,13 ) j^p ^ [ l.flfz), f2(z), 

12^(2), fl(2), f2(z) ]dz 

• -lab 3Uo3v^ 

2 4 c*' <rt> + 


3 u « 9 V r-v 

( t - 9 ^^ , / £i\2 , Aw\2 


+ 2q1 

9 u 

0 

9t " 

l-i . ^ ''o 

3t ■* 'n • 

-Li 

at 


9 u 

3^w a 

2 

+ 2q1 

o 

"Tt “ 

a“tax Tt • 

a w 
atay 

+ It 

/- ^ X'|2 




2 2 
, 3 Wo 3 w <5 

^3tax^ ^3t5y^ 


+ OTq ^ ^ 

3t V am It • 




-V, - M N M N 

(•—) [22 2 2 

m=1 n=1 r=1 s=1 


} } ^ ''mn-Ve»un<’'^- 

^rs-V^ •*up(’') 


^ * °rs- V » + 0„„-2>.„(£) 


mn wm 
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2q1 •’’xs' ’>> 

+ ’ \s‘ ’ *ys*'’* 

^ C^„. «!„„(£) . 0^^ . 1 . 0^(1) 

°nm- ‘ vm^’^ ‘ °rs‘ b ' 

+ \n- ®xm(*) •* xnf' > * \s* l^xrf^ > •»'xsf''> 

+ V- • '^rs- ®yr'^ ’ '’ys'”) 

^2 ^2 • o„„- 0i.<« • Ivn<n) . 0^3. 0-^(1). 

■" • °mn- ®wb'*> •»wn<’>> • “rs* •?«(’n) 

+ » i * ’^mn- »xm(^) •»x„<’’) • °rs- ^’wx^ ^ 

+ E • '^mn. • *yn^’'> * Ops* '*'*’> 

(21,25) 

Now using the equation (2222) and tiie notations for 
integrals given in Appendix-rA, we get a set of linear 
homogeneous, algebraic relations which can be concisely 
put in matrix form as follows? 



CENT^ ''' ’ 


[s3 


-*^ [“]} 


rs = 0 


(2.26) 


This is a standard eigen-value problem, where is the 
non-dimensional frequency parameter. 


^n " 


Ej h 


(2.27) 


]is the stiffness matrix and[|jis the mass matrix the 
elements of which are as follows. 


Stiffness Matrix Elements 


2 11 00 10 01 
a [All. I 1 mr . Jins + 2 A16>^ . Ilmr Jins 

2 00 11 
+ A66. X . Ilmr J1 .ns ] 

2 10 01 11 00 
2 a [ A12. X. I6mr . J6ns + A16 I6mr J6ns 

2 

4 A26.\ . 16®r . J6ns 4 A66.X0 I6mr J6ns 


12 00 


10 02 


= -2a [ I7ffir J7ns C4E11 4 X iTmr J7 bis. C4.E12 


4 C4X E16 (2I7mr 


01 02 
J7ns + 17 njr 


10 

J7ns ) 
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■^3 00 12 o 01 11 

+ X . C4.E26 ( I7nir • J7ns) +X . I7mr • J7ns . C4.E66 ] 


2 11 00 

SI 4 = 2 a [ISmr . J8ns (CIBII - C2 Ell) 

01 10 

+ 2 X ISmr J8ns (C1B16 - C2 El 6) 


2 00 11 

+ X ISmr J8ns (C1 B66 - C2 E66) 3 


2 10 01 

S15 = 2a [X .I9mr J9ns (C1 B12 - C2.E12> 


11 00 

+ I9inr J9ns (C1B16 - G2.E16) 


2 00 11 

+ X I9mr J9ns (C1B26 - C2 E26> 


01 10 

+ Xl9mr J9ns (C1 B66 - C2.E66) ] 


2 2 00 11 01 10 
S22 = a^ [ A22 X I2inr J2ns + 2A26 X I2inr Jrns 



11 

00 






A66 I2inr 

J2ns 3 






02 

10 





= 

-2a [XII Omr 

J1 0ns C4 

E12 





12 

00 


00 

12 


"f 

IlOiar J10 

r% 

ns C4 q^E16 



J1 0ns 

C4.E22 

02 

+ 

X^ C4 E26 f2 I10iar./J10ns + 

11 Omr 

+ JlOns) 


11 

01 





4 

2^. 11 Omr 

JlOns C4* 

E66 J 
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2 01 10 

S24 = 2a [Xinmr J11ns (C1 B12 - C2.E12) 

11 00 

+ I11mr J11ns (C1 . B16 - C2. E16) 

2 11 

+ X I1lmr J11ns (C1 B26 - C2 E26) 


10 01 

+X.I11mr Jllns (C1 B66 - C2. E66) 


2 2 00 1 1 

S25 = 2a [X I12mr J12ns (C1 . B22 - C2. E22) 


1 0 01 

+ 2.X,I12mr J12ns (C1 .826 - C2.E26) 


11 00 

+ 112 ar J12 ns (C1. B66 - C2. E66) ] 

2 2 1 ^ 2 
S33 = a X I3inr J3ns [ A44 - 6D44C4 + 9.F44 C4 ] 


2 

+ 2a .Xj:3n!r J3ns [M5 - 6 D45 . C4+ 9 F45 C4 ] 


2 I • 2 

+ a I3inr J3ns [A55 - 6D55 C4 + 9F55 OT ] 


12 00 2 2 20 02 2 
+ I3iBr J3nsH^^ + 2.X* I3iBr J3ns H12. C4^ 


21 01 2 4 00 22 2 
+ 4 ,Xi3nir J3ns H16 +X , I3iBr J3ns H22. C4 
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+ 13 ^ ^ 2 ^ ^^2 


3 00 01 

S34 = 2 a xI13mr J13ns [Cl A45 - 3 D45 (C2 + C1.C4) 


+ 9F45 Ca. C4 ] 

3 01 00 

+ 2a 113 mr J13ns [A55 Cl - 3D55 (C2 + C1C4> 


+ 9 F55 C2. C4 ] 


12 00 

+ 2a I13ira: J13ns [C2.C4.H11 - C1.C4.F11 ] 


2 10 02 

+ 2a. X . 113 mr J13ns[ (C?.C4.H12 - C1«C4.F12y ] 

11 01 02 10 
+ 2a. A,(2 113 mr J13 ns + I13mr J13n^C2.C4.H16 


- C1 .C4. FI 6) 

3 00 12 

+ 2a X'^ I13mr J13ns (C2e4 H26 - Cl .C4.F26) 

2 01 11 

+ 4a X^ I13nir J13ns {C2.C4.H66-C1.C4.F66) 


3 00 01 

= 2a .X,I14 mr J14ns (C1 .A44 - 3(C2 + C1.C4) D44 

01 00 

+ 9C2.C4 F44) + 2a'^ I14mr J14ns (C1 A45 - 3 {C2+C1C4) 


D45 + 9C2.C4F45) 
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02 10 

+ 2X a I14mr JUns (C2.C4.H12 - F12.C1.C4) 

12 00 

+ 2a 114 nrr J14 ns (02.C4.H16 - F16 C1.C4) 

3 00 12 

+ 2a X I14mr J14ns (C2.C4 H22 - Cl .C4 F22> 

2 / 01 11 10 09 
+ 2a X (C2.C4.H26 - Cl .C4 F26) (2 I14ffir J14ns + I14mr J14ns) 

+ 4aX I14mr J14ns (C2.C4 H66 - Cl .C4 F66) 

4 o 

S44 = a I4mr J4ns (C1 . A55 - 6C1 .C2 D55 + 9C2^ F55> 

2 2 

+ a I4mr J4ns Cl Dll - 2F11. Cl .C2 + C2^ H11 
2 10 01 

+ 2a X I4inr J4ns Cl D16 - 2C1C2 F16 + C2^ H16 


2 2 ' ' 2 

+ a X I4ii!r J4ns (C1 D66 - 2C1C2.F66-i- G2^ H66) 


4 o 

S45 = 2a I15inr J15ns (C1 A45 - Cl .C2.D45 + 9C2^ F45) 

2 2 

+ 2a X I15nir J15ns (C1 D12 - 2C1 .C2 F12 + C2^ H12) 

2 

+ 2a I15mr J15ns (C1 D16 - 2Cie2.F16 + C2 H16) 

2 2 11 2 

+ 2a X I15ffir J15iis (C1 D26 - 2C1C2.F26 + C2^H26) 

2 °1 10 2 o 

+ 2a X IlSmr J15ns (C1 D66 - 2Ct .C2 F66 + C2‘^H66) 



56 


S55 


4 QO 00 
a I5mr JSns (cr 


+ a 




00 


11 


ISmr J5ns 


01 


10 


+ 2a X I5inr J5ns (C1 


11 


00 


+ a I5mr J5ns (C1 


A44 - 6C1 .C2 D44 + 9C2^ F44) 

(C1^.D22 - 2C1 .C2F22 + C2^H22) 
2 

D26 - 2C1.C2oF26 + C2^ H26) 

2 

. D66 - 2C1 .C2 F66+ C2^ H66 ) 


Atoss Matrix El mpn-hc 

00 00 
A^1 1 = P. I1mr , Jins 

M12 = 0 

01 00 

]V113 = 2c^/a . nmr. J7ns 

00 00 
Ml 4 = 2q1 • ISmr • J8ns 

M15 = 0 

00 00 
M22 = P. I2nir , J2ns 

00 01 

M23 - 2q2/b. IlOmr • J1 0ns 


Aii24 = 0 

00 00 
M25 = 2q1*I12mr , J12ns 


00 00 2 00 00 11 
m3 = P.iSffir. J3ns + 12 (i/a I3iar J3ns + 1/b .lamr J3ns) 

01 00 

M34 = 2l3/a • IlSair J13ns 
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Px “) -K (f. 1) = 0 

P y =Py ('.1) = 0 (2.28a) 

These boundary conditions have been used only for crossply 

31 ^ 

plates by Reddy . While for angle ply laminates following 

boundary conditions have been employed* 

Uo,(0,T) ) = u^ (1,T)) = 0 

Vq(S ,0) = Vq U , 1 ) =0 

Wq' ( ^ 0) =. Wo (5, 1) = (0, Ti) = Wq (1^ ) = 0 

P 0) =p^ (fi, 1) = 0 

P y (0,T|) = 0 (2.28b) 

But in the present work we have employed the same boundary 
conditions as given in equation (2.28a) for all types of 
laminate constructions. In Figures 2.3a and 2.3b these 
boundary conditions have been shown schematically. 

2.8 Admissible Functions 

As mentioned in Section 2.6, an admissible function 
should satisfy atleast the geometric boundary conditions 
and should be differentiable half as many times as the 
degree of the system. The boundary conditions specified in 


/ 



59 


the above section can be satisfied by assuming the following 
form for the admissible functions. 


) 

= Cosffi«( S) 

’^un^’^^ = Sin nn( i 


= Sinffiit(g) 

» yyyj('n) = Cosnu in) 


= Sinmn(g) 

» = Sin nic 


= Cosnm(S) 

* ’^xn^’^^ = Sinn^C Tii) 


= Sinm^iCl) 

» Wyj^(‘n) = Cos nuC t) 


Similarly for other boundary condition suitable 

must be chosai* 

admissible functions/ By looking at the expressions for 
the stiffness and mass matrix elements, it can be observed 
that, one will have to obtain corresponding integral values 
for that boundary condition* Thus this method can be 
extended to obtain results for some other boundary conditions 
by using corresponding integrals® 

2 *9 Mode Shapes 

After solving the eigenvalue problem of equation (2®26) 
the deflected shape of the plate corresponding to a parti oalar 
mode of lateral vibration can be obtained by substituting 
the values of the coefficienl^ 0^ in equation (2.23) where 
i indicatos thB node of vibration^ Tha lataral daflectioE 
at any point P(Sj^ is give^ 
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M N 

S 2 

ra=1 n=l 


°mn 


where the admissible functions will take up a form so 
as to satisfy the boundary conditions. Here the 
subscripts k and 1 denote the point in the field. 

Points are selected at intervals of A k and A 1 which 
will depend on the plate aspect ratio a/b. The ratio 
ofA k/Al is always made equal to the plate aspect ratio. 

Using these values of lateral deflection nodal lines can be 
plotted. To get the exact mode shapes for higher modes one 
should consider higher number of terms in the series. 

Moreover depending on the aspect ratio of the plate nuiricer 
of terms in each direction should be varied to get proper 
convergence. 

Thus we have formulated a general solution for the 
problems of free vibration of laminated composite plates. 

This can also be used for the static analysis of laminated 
as well as isotropic plates. This analysis is applicable to 
any combination in-plane and out-of -plane boundary conditions. I 
the next chapter numerical computations have been done 
to show the convergence of the solution, effect of aspect 
ratio, fibre orientation, number of layers, span to thickness 
ratio for isotropic, regular as well as hybrid composites. 

Nodal lines have been ebtained for some representative cases. 
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Fig'2*2a statically equivalent force system on the 
laminates associated with membrane action 



Fig: 2'2b StaticaUy equivalent moment system on the 
laminate due to bending 
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Fig: 2-3 a Schematic representation of the boundary 

conditions convened in the present analysis 





Fig: 2-3 b Sch«iatic rqjresentation of die boundary 
conditions omsidered in Rdl [3fl 
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CHAPTER 

results and discussion 


1 

3 . 1 Numerical Calculatinng 


The present formulation which is developed for the 
analysis of hybrid laminated composite plates has been 
checked by comparing the results for Isotropic and regular 
composite plates for which results are available in liie 
literature. All these computations have been carried out 
on DEO-1090 system using standard routines available for 
numerical integration and Eigenvalue solution. In all 
these cases M = 4 and N = 4 has been considered in the 
summation series of the displacement and rotation functions. 
This will lead to 16 terms in each of the displacement and 
rotation functions thus rendering iiie overall stiffness 
matrix and mass matrix [M] to be of the ordet of (80x80). 
A convergence study has been carried out and it has been 


found that M = N = 4 is the optimal value for Ti*ich 
satisfactory convergence can be achieved. Beyond this the 
storage space required and CPU time becomevery high and from 


the point of convergence it has been found that there is no 

the numerical values obtained# 

Computer programme thus generated is later on used 
for the parametric studies on the behaviour of h^rid laminated 
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composite plates. Three types of materials have been 
considered and non-dimensional mechanical properties 
of these materials areshown in Table 3.1. For the hybrid 
laminate GfRP fibre having higher E|,/E^ ratio has been 
considered since this gives better properties than 
the second type of GFRP for which E^/E^ = 25 . 

For the purposes of tiiis investigation only four 
layers are considered. Numerical results have been 
obtained for a hybrid laminate in which Grar^iite /Epoxy plies 
form the outer layers while Kevlar /Epoxy plies form the 
inner layers to study the effect of various parameters such 
as aspect ratio, length to thickness ratio of the plate, 
fibre orientation and nuntoer of layers. Shear deformation 
effects have been studied by considering thick plates as 
well as thin plates. In what follows we shall discuss 
effect of each of these parameters on natural frequencies. 

3 »2 Convergence Studies 

Convergence study has been carried out for the non- 
dimensional frequency parameter for isotropic, regular 
as well as hybrid laminated composite plates. In Fig. 3.1 
the rate of convergence of the frequency parameter with 
the number of terms considered in the series form of 
the displacement function has been shown, for an antisymmetric 

4 layered angle ply laminate. The rate of convergence 
for the third and -Uie fourth mode have been shown# It 
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can be seen from the Table 3.2 that for M = N =,4 
satisfactory convergence can be achieved even for 
higher modes upto four. This results in an order of 
(8G X 80) for the overall stiffness and mass matrices. 
From Fig. 3.1, it is observed that the nondimensional 
freguency attains asymptotic values "with increase in 
number of terms. 

Table 3.3 gives the percentage change in the 
values of the non-dimensional frequency parameter with 
number of terras in the series. It is clear from this 
table that the change in the value of the frequency 
parameter has reduced appreciablyfor M and N equal to 4 
and that sufficient accuracy has been achieved up to 
the 4th mode. 

3 *3 Isotropic Plates 

The present formulation developed for the free 
vibration, studies of hybrid laminated composite plates 
has been checked by compailing the results obtained with 
those available in the literature. In Table 3.4 a 
comparison of the non-dimensional frequency parameter 
for a simply supported isotropic square plate with length 

j(“lO 

to thickness ratio (a/tt), has been presented with the 
corresponding values obtained by other investigators. 

Since the length to thickness yatid is 10, this is a 
thick plate and hence transverse shear deformation effects 
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can not be neglected* Also for higher modes of vibration 
rotatory inertia effects can be ignored. This can be 
clearly seen by comparing the results obtained with 
those obtained from classical plate theory. The percentage 
error in the classical plate theory increases for higher 
modes. It can also be observed that the results obtained 
using the present solution are in good agreement with 
those obtained from other methods of solution. The 
percentage error when compared with the exact solutions of 
3~D linear elasticity - for the fifth mode is found to 
be less than 1 *2* When compared to those obtained from 
the higher order theory of Kamal and Durvasula , it can 
be seen that there is not much difference between the 
two. This is quite expected since the displacement field in 
both formulations are almost identical. It is observed 
that with as few nuadjer of terms as M = 4 and N = 4 in 
the series form of displacement function convergence Is 
faster than other methods herein more nurfser of terms were 


considered. This is due to the reas(»i that unlike other 
solution methods here, the stress conditions considered 
resemble closely that of i^ree dimensional solution 
method. 


It is well known that transverse shear effects 



reduce the natural freqwem^ejS 
transverse shear 'eff«®i^ 
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has been plotted in Fig. 3.2 for length to thickness 
ratios varying from 10 to 100. It can be seen that 
for ratios ranging from 100 to 30, the transverse shear 
effects are almost negligible since they belong to the 
class of thin plates. Where as after this transverse 
shear effects become predominant. The frequency 
parameter reduceH by 0.6 percent as the ratio is reduced 
from 30 to 20 'vrfiile for 20 to 10 it reduces by 1.86 
percent. From the figure it becomes clear that transverse 
shear effects can not be neglected for plates with 
length to thickness ratio ^ 30 « 

The effect of one more important parameter, that 
is the aspect rati o of the plate also has been studied 
and presented in Fig. 3.3. Non-dimensional fundamental 
frequency parameter for a thin isotropic plate (a/h = 50) 
for aspect ratios (a/b) varying from one to five has been 
obtained. It can be seen toat the frequency parameter 
increases, continuously with Hie aspect ratio. Ibis is 
because, as the aspect ratio increases the stiffness of 
the plate reduces hence the natural frequency goes down. 
Since the non-dimensional parameter has been represented 
as 

the 
that 
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frequency 


parameter values, when they 





are defined as belows 


(3.2) 


where 

X- aspect ratio of the plate. 

Through out this work we have used the non-dimensional 
parameter as defined by equation (3.1). 

3.4 Composite Laminated Plates - Comparative Study 

Table 3.6 shows the values of the non-dimensional 
frequency parameter (K^) along with those obtained by 
other investigators for simply supported antisymmetric 
rectangular angle ply laminated plate. Values of the 
frequency parameter for different modes have been given. 

It can be seen that the values differ for all the modes 
and all the values are slightly higher than those do tained 
by other investigators. This is due to the reason that 
the boundary conditions we have considered are slightly 
different from those considered by other investigators as 
explained in Section 2.7. The particular specifications 
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for the boundary conditions increase the stiffness of 
the plate, ultimately resulting in higher non-dimensional 
frequency parameters. Hence in all these even though 
it is written as simply supported, boundary conditions 
at the particular edges are slightly different. In 
the present work we have considered the boundary conditions 
as given by equation (2*28a)<. 

Effect of Length to Thickne ss Ratio an d Aspect Ratio 

3*5*1 Symmetric Cross Ply Plates 

Figure 3.4 shows the variation of the non-dimensional 
frequency parameter with length to thickness ratio for 
the two types of Graphite/epoxy laminates. For the 
symmetric four layered cross ply plates considered, it was 
found that the non-dimensional parameters were hi(^er for 
the Graj^ite-Epoxy laminate for which the ratio of longitudinal 
to transverse modulus is higher. This is because for this 
laminate the effective stiffness of the plate is higher. 

Here also it was seen that the transverse shear effect was 
predominant as the length to thickness ratio was reduced 

below 2G. 

Figure 3.5 shows the plot of the variation of the 
frequency parameter for regular laminates made of SraFhite/ 
Epoxy plies with higher modulus ratio, Keylar/Epoxy 
plies and hybrid laminate obtained by stacking Graphite/Epoxy 
plies as outer layers and Kevlar/Epoxy plies as inner 
layers. Here also the same type of behaviour with length 
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to thickness ratio is maintained and that the hybrid 
laminate has frequency parameters in between the two 
regular composites. The frequency parameters for the 
Kevlar /Epoxy plies is lesser when compared to that of 
the Graphite/Epoxy plies, while for the hybrid laminate 
the values are closer to that of the Graphite/Epoxy 
plieso 

The effect of aspect ratio on the two types of 
Graphite /Epoxy plates has been investigated and shown 
in Fig* 3.6. The frequency parameter increases with 
aspect ratio and length to thickness ratio indicating that 
transverse shear effects should be considered for thick 
plates. Since the Graphite-Epoxy plies with higher modulus 
gives higher frequencies, we have considered only this 
grade of Graphite /Epoxy to obtain the h^rid laminate* 

For aspect ratios <1, the frequency parameter remains 
constant and in this region the hybrid composite and 
Graphite /Epoxy laminates have almost equal non-dimensional 
frequency parameter values. For aspect ratios > 1, 
the frequency parameter increases linearly thus indicating 
high dependence on aspect ratio# Figure 3*7 shows 

similar behaviour of the frequency parameter with 

I 

aspect ratio for various length to tiiickness ratios for 
regular and hybrid laminates# 
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3»5.2 Antisymmetric Crossolv Laminates 

Figures 3o8 and 3 *11 show the variation of the 
non-dimensional frequency parameter with length to 
thickness ratio and aspect ratio of the plate for anti- 
symmetric regular and hybrid composite cross ply plates* 

In this type of construction plies having fibre orientation 
Qo and 90° are arranged alternatively. The results are 
presented for a four ply antisymmetric cross ply laminate. 

Since the elements of coupling matrices are non-zero, 
their effect is to reduce the frequency parameter values. 

By comparing the results obtained for symmetric and anti- 
symmetric cases, it can be seen that the amount of reduction 
for a square hybrid laminate is 25.3999 for length to thickness 
ratio 10 while it is 29.379^ for length to thickness ratio 20. 
But the behaviour with aspect ratio and length to thickness 
ratio remains same as in the case of a symmetric laminate. 
Figure 3.8 indicates that the frequency parameter reaches 
an asymptotic value for aspect ratios upto 4. For plat^ 
with aspect ratios upto 2 it can be seen that the transverse 
shear effects are almost negligible. 

3.5.3 Antisymmetric Anqleply Laminates 

In this type of laminate construction plies having 
fibre orientation t© and -© are arranged alternatively. 

In this type of laminate construction also, the coupling 
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niatrices B, E and F vill not be zero. However they 

tend to zero as the number of layers is increased. 

Two types of antisymmetric angle plies, 45°/-45°/45°/ 

-45° and 30°/-30°/30°/-30 ° have been considered and 

the variation of the frequency parameter with length 

to thickness ratio and aspect ratio have been presented 

in Figures 3.9, 3.10, 3.12 and 3.13. Even though 

we observe similar behaviour in these cases also, it can 

be seen that the frequency parameter has highest values 

for 45° antisymmetric ply due to the fact that the 

stiffness of a laminate is highest for fibre orientation 

equal to 45°. It can be seen that the frec[uency parameters 

are 6*44 percent higher for a hybrid laminate with length 

to thickness ratio 10 and aspect ratio 1 than a corresponding 

30° antisymmetric laminate. In all Ibese cases it can be 

seen that the hybrid laminates have intermediate values to 

those of the regular composites made of GrajAite/Epoxy 

and Kevlar/Epoxy plies. For thicker plates it can be seen 

that the values for the hybrid laminates are very much 
near to those of the Graphite/Epoxy plies. It can also be 
observed that the transverse shear effects for plates witJi 

same length to thickness ratio increases with aspect ratio. 

In the case of antisymmetric angle plies, it can be seen from 

Fig. 3.9 that the non-dimensional frequency parameter reaches 

an asymptotic value for aspect ratios upto 3 unlike 

anti synme trie crossply laminates where it went up to aspect 

ratios equal to 4« 
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3*6 Effftet of Number of LayA-r.e; 

The effect of number of layers on the non-dimensional 
parameter has been presented in Fig. 3.14 for symmetric 
crossply laminated plate and for a angle ply laminnate in 
Fig. 3*15. It can be seen that for constant length to 
thickness ratio and aspect ratio the frequency parameter value 
increases with the number of layers. But after 12 layers 
it can be seen that the value becomes almost constant. The 
reason for this is that the plate stiffness coefficients 
A16, A26, D16, D26, HI 6 and H26 become negligible with 
increase in the number of layers since, A16 and A26 are 
inversely proportional to the number of layers, D16 and D26 
are inversely proportional to the third power of the number 
of layers, HI 6 and H26 are inversely proportional to the 
seventh power of the number of layers. The non-dimensional 
frequency parameter value goes up by 22.1 89^ for hybrid 
cross ply laminate as the nuiriber of layers is increased 
from 4 to 8, while for 45° symmetric angJe ply laminates 
this value is 24 .l^. 

3 .7 Effect of Fibre Orientation 

Figure 3.16 shows the variation of die fundamental 
non-dimensional freque ncy parameter with fibre orientation 
for square regular as well as h'^rid composite laminated 
plates. From the figure it can be seen that the variation 
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of non-dimensional parameter is almost a normal 
distribution curve with maximum value at 45“ fibre 
orientation since the plate stifness will be maximum 
for this fibre orientation. 



Table 3»1 


Composite Lamina Properties 


Properties 


Composite Materials 

oo 09 1 0 

Material-! Material-!! Material-!!! 


El/Ej 

40.0 

25.0 

14.82 

Glt/E-p 

0.6 

0.5 

0.375 

G»j»7-/Ey 

0.5 

0.2 

0.375 

Glz/Et 

0.5 

0.2 

0.375 

^LT 

0.25 

0.25 

0.34 

Material 

! - Grai^ite/Epoxy (High Modulus 

variety) 


Material !! -Graphite/Epoxy (Low Ibdulus variety) 
Material !!! - Kevlar /Epoxy 



Table 3.9 
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Convergence^^gf the Non-dimensional Freni, 

Number of Terms for Antisvmmetrir ennio D^trSr- ^rameter with 

y ^sniinate + 


M 

N 

3rd Mode 

4th Mode 

Mode 

6th Mode 

7 th 
Mode 

1 

2 

46.5103 

- 




2 

2 

44.3101 

47.6528 

69.1 57g 



2 

3 

43 .541 7 

47.0058 

68.2643 



3 

3 

42.8901 

45.6913 

62o3848* 

62.6739* 

74.1432 

3 

4 

42 .^878 

45.3915 

62.0606* 

, * 






62«4431 

71 .4482 

4 

4 

42 .3430 

45.1760 

60 . 1573 * 

61 .4249 

70 .3716 
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Aspect ratio (a/b) 


K =(,) 
n n 




ic 


,b2(P/Ejh2f 


1 

2 

3 

4 

5 


5.96429 
14.9164 
29 .8739 
50 .7440 
77.6303 


5.96429 

3.7291 

3.3193 

3.1715 

3.1052 



T5 

d) 

U 

<D 

(X5 

u 

D 

O 

tl4 

fO 

o 


vD 

o 

coi 

0)1 

HI 

<0, 

H 




CM 


fO 


li 


0 

-P 

O 

S 

0 

P 

0 

Oh 

>- 

o 

c 

0 

D 

cr 

0 + 

U 0 
Ll, ^ 
0 

r ^ ,--1 

0 Ol 
C 


o 

•H 

<f> 

c 

0 o 

B a 


0 

+5 

•H 

0 


♦H 

T3 


^ o 
c 


o 

2: 

O 

d 

o 


0 


TJ 

0 

-p 

p 

O 

a 

a 


o H 

O if} 


T3 

vO 

N 

CM 

N 

Kf 

CD 

lO 

CM 

C 

a 

« 

a 


0 

00 


O 


+> 

P 

0 

CQ 


CO 

ID 

lO 


\ ^ 

•p o 
C 

0 -P ^ 
0 23 X 
0 H ^ 
P Ow 
Oh 0 


*0 

d 

0 


r- 

CM 


0 


H 
0 
s 

0 _ 
Q 


CM 

CM 




U 

o 

0 

jd 

CO 

2: 

d 

o 

T3 

0 

0 


N- 


0 

tn 


CM 

X 

CM 


0 

te. 


O du-p 


c 

o 


•p 

0 

d 

C3> 

•H 


VO 


Ov 

o 

vO 

CO 

O 

o 

CO 

r— 

IN 

tO 

00 


N 

lO 

CM 

vO 

CO 

CO 

T— 

CO 

• 

• 

• 

•' 

• 

CO 

CM 

o 

o 


CM 



vO 

N 


D 

vD 

T~* 

o 

o 

t— 

0 

O 

CM 

vO 


o 

0 

O 

to 

lO 

lO 

CM 

> 

• 

« 

« 

• 

• 

P 

a\ 

to 

<r- 

vD 

to 

D 


CO 

to 

lO 

N 


CO 

to 

o 

N 

lO 

o 

ON 

CO 




NO 


« 

1 ♦ 



to 

lO 

N 


CO 

lO 

NO 


CD 

• 

CM 

00 


00 

CM 

CO 

« 

00 


vO 

IN 

• 

CO 

O''- 


lO 

CM 

N 

• 

00 


0 

OQ 

0 

•p in 

ON 

lO 

N 

to 

N 

'*W' 

0 11 

lO 

CO 

NO 

, Hr- 


s 


CM 

9 

lO 

• 

1 ^ 

1 0 

C^> 

• 

to 

9 

tn 

c 

00 

to 


o 

o 

tn 

<N 


CO 

to 

N 

o 


23"'' 

o 

CO 



t- 

to 

€) 

to 

•H IN 

to 

IN 

▼- 

w 

>o 

N 


0 0 p 

0 -P C 

• 

to 

« 

CO 

9 

• 

00 

9 

IN 

• 

.r-1' 

a 

rS' 

0 0 0 
■ H H , x: 

CM 

to 

Ov 

On 


CM 


04 CM lO ^ 


80 


0 

H 

-P 

u 

0 

a 

o 

u 

a 

H 

0 

•H 

P 

0 

-P 

0 


P 

o 

%4 

CO 

(• 

CO 

0 



81 





Fig; 3>2 Variattofi of non dimension^ frequency parameter with lengtt) to 
thickness ratio for«i isotropic s<|uare i^ate. 



« .S;' 










Length to thickness ratio , ag1» 

Variation of non (limenfti<»iai #re<|ijency parameter with tength to thickness 
ratio for ^mmetric cross ply s<^are plates. 





Non diinensionat frequency paramettr k 
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Non dimensional frequency parameter 





Rg: 3-8 fret^iency pwsffindQ^ v4ih length-b ihickTess ratio 

for ariisywiiiiliiic cross ply laminated plates 
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Variation of non-<&nensioncft frequency parameter with 


aspect ratio for antisymmetric cross pty plates(Gjg(M/9d 
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Fig: 342 Variation of norMlimensionM. frec|uency paramat^ with 
aspect ratio for antisymmehic plates 
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Pig: 3*14 Variation of non-dimensional frequency parameter with ixmber 
of layers for symmetric cross-ply plates 

> ' . ' > i' ; V ■ 
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CHAPTER 4 



» . 

Use of hybrid laminated plates is an effective 
means to achieve intermediate properties and to tailor 
properties as per design requirement using various para- 
meters such as stacking sequence, nun4)er of layers, fibre 
orientation etc. From llie limited studies carried out 
on the free vibration response of regular as well as 
hybrid laminated composite plates, it can be concluded 
that 


(1) Rayleigh-Ritz energy methdd can be successfully 
employed for the free vibration analysis of regular 
as well as hybrid composite plates, ^preciable 
convergence of liie solution can be achieved with 
as few as 4 terms in -toe summation series thus 
reducing the computational time. 

( 2 ) Tlie properties of the hybrid laminate are intermediate 

% 

to those of single material coiq>osites. 


Cd) The transv^se shear effects can not be neglected for 



plate® length to thickness ratio less than or 

'■ Kt- .4. 

eq»*l to 120 and rotatory inertia effects can not be 
■neglected for higher modes. 

'I"'*'--.’' 

iTTile'hataral frequency of vibration reduces as the 
plate aspect ratio increases. 
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(s) For laminates with any stacking sequence fibre 
orientations + 45° give the highest fundamental 
frequency# 

Scope for further extension 

The following is the list of several possible 
interesting aspects for extending the present work 

(1) For free vibration analysis different boundary conditions 
can be considered. It is very important to investigate the 

free vibration response of plates with mixed boundary 
conditions* 

(2) Transverse shear and rotatory inertia effects on 
trapezoij^al as well as skew plates can be studied. 

(3) The effect of delamination on natural frequencies can 
be investigated., 

(4) The present formulation can be used to study the 
buckling and bending response of hybrid composite 
plates under various boundary and loading 


conditions 
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appendix a 


The list of the integrals involved in this formulation 
are as follows: 
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For the simply supported conditions specified in Section 2*7 
and admissible functions given in Section 2.Sfthe values of 
the integrals have been listed. These integrals have been 
obtained by numerical integration technique using available 
routines in the DEC-1090 system. 
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